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0^ ' six (including both gauge and fermion operators) which must be added to the unimproved action 

to form an 0{a^) improved action. Among them, the total number of dimension-6 gauge operators 
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I. INTRODUCTION 

In recent years, there has been a surge in developing and applying improved actions for the numerical simulations of 
lattice QCD. Up to now, most research has focused on improvement of Wilson fermions in an effort to reduce the OJa) 
cutoff effects in the simulations. On the other hand, the absence of 0{a) errors for the staggered fermion action pp] 
and the complexity of the staggered formalism mean that its 0{a^) improvement has received little attention. Almost 
ten years ago, Naik proposed adding a third-nearest-neighbor term to the standard staggered fermion action to 
remove some O(a^) effects ^. His study was based on the Dirac-Kahler equation, not on the standard staggered 
formulation. Although these two fermion formulations are the same in the free case, they are quite different when the 
gauge interactions are included and the difference is of the order of a^. So, the Naik term may not remove all 0{a^) 
errors from the simulations by using staggered fermions. This statement was demonstrated by the recent numerical 
simulation from the MILC group [Qj^] 

One approach to improvement is to construct a perfect action |^. The classical perfect action for free staggered 
fermions was already proposed in ref. M. In this paper, we will apply Symanzik's improvement scheme to staggered 
fermions and discuss its on-shell improvement through 0{a^). We will show that including only the Naik term in the 
improved staggered fermion action is not enough to remove all order a^ errors from on-shell quantities. Meanwhile, 
because both the standard Wilson gauge action and the standard staggered fermion action have 0(0"^) errors, we 
must improve both of them at the same time. We will show that these two improvements are not independent but 
connected by an isospectral transformation of the gauge fields. The recent calculations of the MILC and Bielefeld g] 
groups can be easily explained by the result of our analysis. 

This paper is organized as follows. In section R we will discuss the 0{a^) improvement of the stagge red fermion 
action by finding all linearly independent dimension-6 operators following Symanzik's scheme. In section |lll| , we will 
give the coefficients in tree level up to order a^. In section fS/\ , we will give the general isospectral transformation 
of fermion fields and gauge fields and give the form of the simplified on-shell improved action. Section ^ is the 
conclusion. The computations of the tree-level coefficients are presented by two appendices. 

II. POSSIBLE COUNTERTERMS 

When Symanzik's improvement scheme is applied to construct an 0{a^) improved lattice action, the first step is 
to find all dimension six operators which are scalars under the lattice symmetry group. These operators, treated as 
counterterms, are then added to the action to remove all 0{a^) errors from physical quantities. Before doing that, 
we will introduce some notation which will simplify our presentation. For the transformation properties of staggered 
fermion fields, the reader is requested to consult ref. M, which we will refer to as paper (I) in the following. 

A. Definitions and Notations 

The whole lattice can be viewed as being composed of elementary hypercubes consisting of 16 lattice sites. We will 
use X to label the individual lattice sites and y, which has only even coordinates, to label each hypercube. A site 
inside a hypercube is represented by a "hypercubic vector" A, whose components can only take the values of either 
or 1. The relationship between these 3 vectors is given by 

x = y + A. (1) 

The hypercubic fields are defined as: 

XA{y) = ^xiy + A-), (2a) 

XA{y)^lxiy + A). (2b) 



With the notation 



xMill)x=Y. XA{y)MiU)^sXB{y), (2c) 

y,AB 



and 



(7S «) S.f)ab = ^Tr(7;^7s7s7F), 
we can write the standard staggered fermion action in a compact form as 



Sf = X 



J2 (7p ^ n'D^, 



m 



Furthermore, when discussing the fermion operators, we wiU use the foUowing notation: 



(2d) 



(2e) 



(3a) 



^' = E(^a<®^)^a<(7m®^)^m: 



(3b) 



^' = II (7m ® ^)^M (7/. ® ^)2?M (7m ® ^)2?M, 



(3c) 



P:^!? = ^ (7^ ® /)I?^, (7, ® /)!?, (7^ ® /)I?^ 



(3d) 



B. Fermion bilinears 



The lattice symmetry group of staggered fermion action 19U2| includes: translation, reflection, rotation, charge 
conjugate and a continuous C/y(l)- When the mass parameter m is 0, there is a second continuous Ua{^) symmetry. 

When 771 = 0, we can identify the following five independent operators which are scalars under all symmetry 
transformations including ^7^(1): 



Oi = x2?'x, 



(4a) 



02^x^{v^V-V'D^)x, 



(4b) 



03^x^{v^V+V'D^-2V^)x, 



O^^X {V^Tp+TpV^ - 2V'^V) X, 



(4c) 
(4d) 



05=xV>\- (4e) 

For the case of non-zero fermion mass, the ^/^(l) symmetry is violated, and there are two more allowed counterterms: 



06=mxV X, 



(5a) 



Oy = mx'D^X- 



(5b) 



C. Four fermion operators 

When considering 0{a^) corrections to the fermion action, we must examine not only dimension-6 operators bihnear 
in the fermion fields, but also four fermion operators of dimension of six. 

Using hypercubic coordinates, we can connect the spin and flavor indices in staggered four fermion operators in 
combinations of the form: 

{xMiU)x)' = J2T.^^(y^^(^^ABXB{y) Y.xc{y)M{U)cj,XD{y). (6) 

y AB CD 

However, the color indices in such an operator might be combined in four ways: 

XaSaa' (m (W)x) ^, Xbkb' (m{U)x) ^, , (7a) 

XaSab' (M{U)x)^^ XbSba' (M(Z^)x)^,, (7b) 

Xatl,a' {m {U)x) ^, Xbtl,, (m {U)x) ^, , (7c) 

Xai^ab\M{U)x)^Xbtla\M{U)x)^^, (7d) 

with 

t^ = ^, (8) 

where A' are the SU{2>) Gell-Mann matrices, and as usual, the repetition of the indices a,a',b,b' and i means sum- 
mation. Because of the completeness relation of the matrices A* 

8 . . 2 

/ , Ka'Kb' = 2(5afc'(56a' " -zSaa'Sbb', (9) 

1=1 

the operators with the form of Eq. (Q) and Eq. (^) can be expressed as linear combinations of operators with the 



form of Eq. ( |7c[ ) and Eq. (|7d|). Furthermore, the operators with the form of Eq. (7d) can be expressed in terms of the 
operators with the form of Eq. dTq) by making a Fierz transformation. Hence, we need only consider the operators 
with the form of Eq. (0).[| 

For convenience, we will not write out the links explicitly in the remaining part of this section unless there would 
otherwise be confusion. After applying the staggered fermion symmetry transformation including rotation, reflection, 
charge conjugate and the continuous C^y(l) x [/^(l), we found there are 18 operators which are invariant: 

Ti - (xt'IT^x) ' - (xi"(75®e5)x) ' + Y. [{xt'hi^^I^x) ' - (xf'J^U^I^x) '] , (10a) 

At 

T2 = (xt°Cr®T)x) ' - (xt"(75®C5)x) ' - E [(xi"(7A.®eM)x) ' - (x^"(75M®e5M)x) '] , (10b) 

•^3 = 5I(^^''(^5m®/)x)', (10c) 



*If additional, explicit flavors of staggered femiions are introduced, Fierz symmetry cannot be used a second time, so we 
will need to introduce both flavor adjoint and singlet fermion bilinears, effectively doubling the number of flavor singlet, four 
fermion operators that must be considered. 



•^4-^fo"(7M<»^5)x) , 



(lOd) 






(lOe) 






(lOf) 






(lOg 



^8 = 5I(x^"(75m<»^M-)X 



•^9 = E(^^"(^®^5m)X 



(lOh) 



(lOi) 



•^io = 5l(xi"(75<^eM)x) , 



(lOj) 



^ii-E(^*"(>®^) 
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X) , 



(10k) 



^12= Y. (xi"(7M-<»e5A)x) , 

M#!^#A 



•^13-E(^^"(^®^^ 



(101) 



(10m) 



•^14= ^ (xi"(75M'»e.A)x) , 



Mt^I^T^A 



(lOn) 



•^i6 = 5](xi"(7/..'^C..)x)', 






(lOo) 



(lOp) 



(lOq) 



•^i8 = E(^^"(^M®C^-)x) • (lOr) 

These operators are not invariant under translation. However, the additional terms generated by translations are of 
0{a). For any operator listed above, we can combine it with some higher dimension operator so that the new operator 
is invariant under translation. This new operator differs with the old one by an 0{a) term and both of them have the 
same continuum form. Therefore the translation symmetry does not reduce the number of invariant operators here. 
After adding some higher dimensional terms, we can make the 18 four fcrmion operators listed above invariant 
under translation and rewrite them in terms of the fields xi^) s-i^d x(a;). First, 

•^1 = E X(a^)i"x(a^) E X(a^+e)i"x(a^+e), (11) 

x,a e 

where the sum over e is a sum over the 8 possible lattice displacements of length "1". Second, 

^2^J2 X{x)t''x{^) J2 X{x + v)t''x{^+v), (12) 

a;, a v 

where the sum over v is over the 32 possible lattice displacements of length "\/3 " . Next, 

•^« = EEc;(^)^E^(^)'?5(^)^m^^wcm(^+'^) (13) 

x,a fi c 

i = 3,---,10. 

This equation contains a number of new elements which we will now define. The sum over c is a sum over the 81 
displacements with coordinates c^ = — 1,0, 1. The weight is: 

4 

w{c)^l[{2^\c,\). (14) 



The fermion bilinear C2(x) is given by: 



C^(x)=x(:^)t"Ex(^+«), (15) 



where the sum is over the 8 possible lattice displacements of length "v3 " which are perpendicular to /t direction. 
The phase factors P/i (c) are defined by 

^f(c)-^M(c), Pl''Hc)=eic)rj,ic), 

Pf (c) - e{c)T^{c)v^.{c), Pl^\c) = T^{c}r^^{c), 

PPic) ^ r,{c)Uc), P^'\c) = e{c)r,{c)Uc), 
Pf(c)=e(c)C^(c), p(i°)(c)=Cp(c), (16) 



where 



^A<(^) = ^E(-i)"^ 



■ tJ.\"/ 


3^v 


^/ I 


V^{c) = 


(_l)c,+- + 


1 


Uc) = 


(_1)Cm+i + - 


••+C4 
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e{c) 


= (-1)^^+- 


••+C4 


VI 


4 

>(c) = ll^ 


7m(c)- 



(17) 
^=1 



The remaining 10 operators can be written as: 



i = 11,---,18, 



256 



with the fcrinion bihnear 



B;{x) = I x{x)t''x{x+fi) + x{x)t''x{x-fi) (19) 



and the phase factors are given by 

Pl^'\c) = e(c)C,(c), Pj^^\c) = e(c)r,(c)C^(c), 

We have now discussed all diniension-6 fermion operators which are invariant under the lattice symmetry group. 
Therefore, we can write down a suitable 0{a?) improved staggered fermion action as: 

7 18 

SF=x{V + m)x + o?Y. b,{glma)0, + a' ^ ^^(^o, ma)J^^■ (20) 

i=l i=l 

The reality of the action requires that 62 is imaginary, and all other b' and 6's are real. 

D. Gauge fields 

The Symanzik improvement of the gauge theory action was studied more than a decade ago |lfl,0] . It was found 
that there are three independent six-link products which must be added to the original Wilson action to form an 
0{a^) improved gauge action. The improved gauge action can be written as: 

3 
SG[U]=J2c.{gl)C, (21) 

where the link product Ci are defined as: 

-^0 = ^^ I^eTr(l- Q), (22a) 



/3' 



A = ^El^^Tr(l-^^), (22b) 



^2 = fE RcTr(l-f.i..i), (22c) 




A = fEl^eTr(l- JII 



(22d) 



where ( ) implies an average over orientations. The four c^'s satisfy the normalization condition: 

Co(5o') + 8ci(5o') + 8c2(.9o') + 16c3(5o) = 1- (23) 

For on-shell improved pure gauge theory, it was shown that we can set €3(50) to zero by a change of field variable 
in the path integral. However, we have to be careful when we discuss an improved action which includes the quarks, 
because the change of gauge field variable will also have an impact on the fermion action. We will discuss this issue 
in the latter part of this paper when we discuss the isospectral transformations. 



III. TREE LEVEL IMPROVEMENT 



A natural way to do the tree level improvement is to expand the lattice action to order a} and to adjust the 
coefBcients hi so that the difference from the continuum Lagrangian is of order of a?" . This also improves the free 
fermion propagator through order of a^ . 

Define the gauge covariant hypercubic fermion fields as 



^Pa{v) ^Ua{v)xa{v), 



(24a) 



<^A(y) ^XA(.y)U\{y), 



(24b) 



where lAA{y) is the average of link products along the shortest paths from y to y + A. For the classical continuum 
limit of the standard staggered fermion action (see Appendix |AJ), we find: 



•'V AB *- M 

+a Y^ \igo Y ^aIt^T^TJ^^Fa^ - (75 ® ^5^,)ab^, 



-^E(^m®^)as^: 



-^ffoa' \Y ^aA,(7^ (» T)^g [D„ Fa^] 



- Y M15^^5^)ab {[Df., F^^] + SFx^D^)] Lb (y) 
fj,\ ^ 

where 13^ is the continuum covariant derivative, and F^,^ is the continuum field strength. 
Now, let us define a new set of fermion field variables 



expf-a^^A-DAJ ipA, 



(25) 



(26a) 



Da = y^A exp 



(-E 



AxDx 



(26b) 



where D is defined in Eq. (A3a). This definition is an obvious extension to higher order in a of the order a transfor- 
mation discussed in paper (I) chosen to remove extraneous terms of 0{a) which appear when the original staggered 
fermion action is written in terms of hypercubic variables. The continuum limit of the staggered fermion action in 
terms of the ch field can be written as 



p 2 

Sf= Y My){Y (T/^ ® I)AB[Df. + ^Dl] +mJTW7)AB}My)+oia^). 

•'y AB u 



(27) 



From this equation, it is easy to get the tree- level values of the coefficients occurring in Eq. (^0|). We obtain: 

&i(0,ma) = -i, (28) 

and all other b and b' are zero. 



The coefficients in Eq. (21) were given in references |10| and |ll|. Their values are: 



CO(0) = | Ci(0) = -^, C2(0)=C3(0)=0. 



(29) 



The tree-level values of the four-fermion operators are also calculated (see Appendix ^|^ They are: 



b'i2 = ^, (30a) 



bU = f|, (30b) 



b'ls = f|, (30c) 



and all other 6,' are zero. 



IV. ON-SHELL IMPROVEMENT 

The on-shell improved action is not unique. Given one improved action, we can obtain another one by a trans- 
formation of the fields. However, all these actions are equivalent because all such actions will give the same value 
for a specific on-shell quantity. Thus, we can choose to minimize the number of operators occurring in the on-shell 
improved action by an appropriate definition of the field variables. 

A. Isospectral transformation on fermion fields 

To simplify the improved fermion action Eq. (EQ), we consider the following transformation: 

X — ^ (l + a^eim-^ + a^SiTp^ -f a^ezV^^ x, (31a) 



X — 'X 



(1+ a^e[m^ +0^6^^^ + a^s'^v'^]. (31b) 



This is the most general transformation of the fermion fields x and x which preserves their transformation properties 
under the lattice symmetries. After rescaling the fermion fields and redefining the fermion mass parameter 

m' ^m{l-a'^m^{ei-e[)), (32) 

to the first order in the e^'s and second order in a, we end up with the change of the action 



6Sf — a 



(4 - £3)02 + (4 + £3)^3 + (4 + £2 + 4 + £3)05 

+ (ei - 4 +62 + 4)^6 + (£3 + 4)07] + Oia^). (33) 



The reality of the transformed action requires that 

4 = 4, (34a) 

£2 + £2 = I'^al, (34b) 

and 

El — e[ = real. (34c) 



^We thank G. P. Lepage for pointing out the existence of tree-level contributions to these terms. 



Thus, we can always choose appropriate values of the e^'s and e^'s to make the coefScients &2, ^13, 65, feg in Eq- ( pOJ ) 
vanish after it is rewritten by the new field variables. For example, 

4 = -^(&3 + 62), (35a) 

£3 = -^(&3-62), (35b) 

£2+62 = 63-65, (35c) 

El - e'l = 65 - &3 - 6g. (35d) 

Notice that this argument is valid to any order of g^. Hence, to all orders in perturbation theory, we can always 
choose 

b2{gl,ma) = bsig^, ma) = 65(50, "'«) == higa^ma) = 0. (36) 



B. Isospectral transformation on gauge fields 

The general form of the gauge field transformation which changes the action at 0(0^) can be written as: 

C/,,(x) — . U'^ix) = exp[eX^(a;)] C/^(x). (37) 

U' [x) and U^{x) must have the same transformation properties under lattice symmetry group. The general form 
of X^{x) has been given in ref. ||ll| for the case in which ^^(.t) depends on only the gauge variables. It is the 
anti-hermitian traceless part of another field Y^ (x) : 

X^{x) = Y,{x) - Y,{xy - j^Tr[Y^{x) - r^(x)t], (38) 

with N = 3 for SU{3) gauge theory, and 

-U^ix)Ulix -v + (i)Ul{x - v)U,{x - v)) . (39) 

Under the field transformation Eq. (|37|), the path integral is invariant 

"[df/]exp^-[^l+^«[^l = /'[df/']exp^-[^'l+^'^[^'l 

= f[dU] (1 + AJ + e[ASG + ASf]) exp^^I^l+^^I^l +0{e^). (40) 

The change of the fermion action in Eq. ( pO[ ) is: 

ASp^^Oi + Oia'). (41) 

The jacobian 1 + A J and the change of the gauge action ASq will not generate new terms, but only change the 
coefficients Ci in Sg[U]. Liischer and Weisz |1^ discussed the changes of these coefficients in detail. Their results 
showed that we can choose an appropriate value of e to set 

C3(5o) = (42) 



to all orders in perturbation theory. From Eq. (pl|), we see that if we want this to persist, the coefhcient 64(50: "^^) 
will not be zero in general. On the other hand, we can also set bi{gQ,ma) = 0, however, C3((7o) ^^^^ then in general 
be non-zero. 

In order to find out the possible redundant four fermion operators, we will generalize the argument given in rcf [[l2| . 
The lattice action is rewritten in a concise form as: 



s^J2^xhAU)x + TTMu)] 



(43) 



Consider the small change of a link variable 

5Uf,{x) = ieap{x + c) Re [t''Uf,{x + c)]af3 t" Uf,{x), 



(44a) 



SUUx) = -ze„/3(x + c)Re[i"C/^(a: + c)]„/3 f/^(x)i". 



(44b) 



where e is a real small number. The jacobian differs from 1 in order e^. To the first order of e, we get the following 
identity 

dU^{x) e^l'^'^'^1 lxt''U,,{x)^^^^^^^^xSo.t3{x + c)RerU^ix+c)U 



dU^{x) 



-Tr 



eU,,{x) 



dUU^ix)) 



dU^{x) 



eap{x + c) Re [t°-U^{x + c)]ai3 



Similarly, we get 



dU^{x) e^I^'^.t/] \^ta^^^^)^h^Ml^,^^^^^^)i^^t-u^{x + c)U 



-Tr 



t"f7p(x) 



dIAU^.{x)) 



dU^{x) 



Sapix + c) Im [t''U^{x + c)]ap 



Combining the above two equations, we get 



s[x.xm j ^,arr (^^^!h^^EjMl^e^^(^x+c)[eu,{x+c)U 



xt'^u.Xx)- 



dU^ix) 



+Tr 



'''"' dU,{x) 



f'U^ix) 



eap{x+c) [fUf.ix+c^p } = 



(45a) 



(45b) 



(46) 



In the above equation, we replace eap{x+c) by a x"{x + c)x (x + c+fi), multiply it by a combined phase factor 



18 



Y, <V,i^)PJ.ic), 



(47) 



and sum it on the hypercubic vector c. After substituting hx{U) and fx(U) by the actual staggered fermion action 
and the gauge action, we see that we can add the following terms 



18 



A5=^^e:.F. + aM,04 + 0(a3) 



i=ll 



to the action without changing the path integral to the order of 0{a'^). Notice that because of the identity 



E'?M(c)i^;(c) = ^Mi> 



i = 11,---,18 



(48) 



(49) 



there is only one dimension-6 bilinear in Eq. (Eq). Therefore, we conclude that the four fermion operators Tn ■ ■ ■ Tig, 
whose fermion bilinears consist of two sites separated by one link, are redundant and their coupling constants b'^ ■ ■ ■ b'lg^ 
can be set to be zero. Again, the coefficient of operator O4 will get a change accordingly. Is it possible to get rid of 
some other four fermion operators? From the above discussion, we conclude "no" . Because in the original staggered 
fermion action, x s-nd x ^'^'^ separated by one link, there is no way to generate an operator like J^i, • • • , J-'io when we 
multiply that link variable by some Grassmann variables. 



10 



C. 0{a^) on-shell improved action 

The 0{a^) improved action for lattice QCD can be written as: 

Sqcd^Sg[U] + Sf[x,X,U], (50) 



where Sp is given by Eq. (gO|) and Sq is given by Eq. (|21|). 

The operator O4 involves the lattice sites in two nearest neighbor hypercubes. There exists another operator O'^ 
which only involves lattice sites inside one hypercube and is equivalent to O4 up to an operator of dimension of seven. 
We can construct this operator by replacing the link variable Uf^{x) in the Dirac operator "p with a modified link 
(e.g. the MILC "fat link") W^,{x)U^{x) §]: 



01 = E X(^)^ [W,{x)U,{x)x{x + A) - Ulix - fi)Wlix - fi)xix ~ A)] • (51) 



This new operator 0\ obeys all the staggered fermion lattice symmetries. The factor yV^(x) has the form: 



+Ul{x - y)U^,{x -i))U,{x + fi~ OjUUx) - 2] , (52) 



and its continumn limit is igQa[D^, F,^^]. Thus, O'^ has the same continuum limit as O4. Because it is simpler, we 
will choose O'^ instead of O4 to make our on-shell improved action. 

From the above discussion, we can see that one possible on-shell improved action for lattice QCD can be constructed 
as 

2 

^QCD — /_^^ii9a)^i + ^F ' (53) 



i=0 



with 



Sp^ =X(P + "^)x + a'^bi{g'^,ma)Oi + a^b^igl, ma)0^ + a^bj{gl,'ma)0'j 

10 
WY,b'M,ma)T,. (54) 



Another possible choice would be: 

3 

^QCD — Z_^ '^i\9o)^i + Sp , (55) 

i=0 

with 

10 
4'^ =X¥ + m)x + a'bi{glma)Oi + a^rigl ma)07 + «' E Kiiglma)^,. (56) 

i=l 

In either case, the tree- level improved action is the same because both 64(0,7710) and 03(0) vanish. 

D. Formulae arranged for lattice computation 



For the operators appearing in Eq. ( |53| ) or Eq. (g5|), the gauge part, given by Eq. (|2l|), is already in a form handy 

for numerical simulation. However, the fermion bilinears in Sp or Sp are represented in terms of the hypercubic 
fields. In order to do numerical calculation, it is convenient to rewrite them in terms of the original field variables x 
and X- 



11 



Let us first consider the four fermion operators. These 10 operators, described in section || and IV, are aheady 
written in terms of the original variables x ^nd %. However, they can be simulated using known Monte Carlo techniques 
only if auxiliary Yukawa fields are introduced so these four fermion operators can be written in a bilinear form. The 
resulting gauge- Yukawa fermion action will be quite complicated with hermiticity properties that depend on the sign 
of the original four fermion coefficients. Further, the positive-definite character of the staggered fermion action may 
be lost unless these new Yukawa terms possess some additional (staggered fermion) C/a(1) symmetry. Since these four 
fermion terms are not present in tree level (or tadpole improved tree-level) approximation, we will not consider this 
question further in this paper. 

Only considering the fermion bilinears, the fermion action Sp in Eq. (p4) can be rewritten as: 



with 



S^P^a'Y^xi^) 



X,fl 



2a 



U^{x)x{x + A) - Ul{x - fi)x{x - /i) 
+ (l - a3{gl,ma))m a^ ^ xix)x{x) 

X 

+a3(.9o: ma)m a^ ^ xix)^^ U{x, x + 2fi)x{x + 2/i) + U{x, x - 2/i)x(x - 2/i) 



-ai(go,TOa) a'^^x(a 



a:, /J, 



6a 



U{x, X -I- 3/i)x(a; + 3/i) — U(x^ x — 3fl)x{x — 3/t) 



Uf^{x) = (l + ai{gl,ma) - 02(90, ma))Uf_,{x) + (^2i9o,ma)Uf^{x), 



(57) 



(58a) 



and 



U^{x) = g XI [u,^ix)U^,ix + v)Ul{x + ji) + Ul{x^ v)U^,{x - i))U^ix - v + ji) 



f^p, 



U{x, x + 2fi)^ Uf,{x)Uf,{x + fl), 



(58b) 



(59a) 



U{x,x-2fi)^Ul{x-fL)Ul{x^2il), 



U{x,x + 'ijl) = Uf,{x)Uf,{x + fi)Uf,{x + 2fi), 



(59b) 
(59c) 



U{x, x-m^ Ulix ~ fiplix - 2fi)Ul{x - 3A). (59d) 

The three parameters ai, a2 and as are real numbers and related to the parameters &i, 64 and 67 of Eq. (|5J) by: 

3. 



ai = -j&i, 



(60a) 



^2 = ^4, 



(60b) 



"3 = 2^7, 



(60c) 



:(2) 



The form of Sp , when rewritten using the fermion fields x ^nd x, is the same as Eq. (pTh but with 02(^0) "^'^) — 



^(2) ^ _5.(1)| 



(61) 



At the tree-level, we have 
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ai(0,ma) = -, 0:2(0, ma) = 0:3(0, ma) = 0. (62) 

8 

One way to do a better job than tree- level improvement may be to use tadpole improvement. Following the work of 
Lepage and Mackenzie |l^, we can replace all links U^{x) which appear in lattice operators by —U^{x) (Of course, 
the normalization of the coupling constant will be different). The constant uq is the mean value of the link matrix 

and is measured in the simulation by the quantity mq = [Re (iTrf/n)] * . The tadpole impr oved tree level coefficients 
are the same as above except that the parameter oi which appears in Eq. ( |57| ) and Eq. (58a) must be replaced by 
two different paramet ers, q° and a\. The value of the parameter a° appearing in Eq. (|57) is Oi = l/Swg while the 



parameter a\ in Eq. (58a) is still 1 



V. CONCLUSION 

Using Symanzik's improvement program, we have discussed the 0{a?^ on-shell improvement of the staggered fermion 
action in a systematic way. Our first step was to find all dimension-6 lattice operators which are scalars under the 
lattice symmetry group. We found that there are 5 linearly independent fermion bilinears that are invariant under all 
lattice symmetry transformations. When the mass parameter is not zero, the ^7^(1) symmetry of the staggered fermion 
is violated, and there are 2 more fermion bilinears that violate only this Ua{^) symmetry and are proportional to 
the mass of the fermions. For staggered fermions, we observed that there are 18 independent four fermion operators. 
Therefore, we have at most 25 fermion operators which can be added as counterterms to the standard staggered 
fermion action to remove all 0{a?') errors from all physical quantities. Including the 3 independent dimension-6 gauge 
operators, we end up with 28 counterterms for the 0{a?') improved lattice QCD with staggered fermions. 

For on-shell improvement, we can use the isospectral transformation of the field variables to eliminate all possible 
redundant operators. Including such field transformations, we concluded that we need at most 15 independent lattice 
operators of dimension-6 to construct the 0{a?) on-shell improved lattice QCD with staggered fermions. Ten of these 
are four-fermion operators, which are absent at the tree-level, and hence of the order of 0{gQa^) at most. The other 
5 are fermion bilinears and gauge operators and only two of them are nonzero at tree-level. Two possible improved 
actions are given by Eq. (|5|) and Eq. (^5|). 

Thus, we found that the Naik term is not the only term in the improved staggered fermion action. It is worth 
emphasizing that to remove the 0(0^) errors from lattice computation, we must use both an improved gauge action 
and an improved fermion action at the same time, not just one of them. 

The recent numerical results from the MILC |0,p| and Bielefeld S groups are consistent with our analysis. Further- 
more, in the free case, our result is the same as Naik's, as should be expected given the equivalence of free Dirac-Kahler 
and staggered fermions. 
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APPENDIX A: CLASSICAL CONTINUUM LIMIT OF STAGGERED FERMION ACTION 

In this appendix, we will evaluate the continuum limit of the staggered fermion action, 

SF = a'^Yl X{x) [^ rif,{x)—{u^{x)x{x + jj) - Ul{x ~ fi)x{x - fi)j + mxix) , (Al) 



■'■Thanks Attila Mihaly for pointing out this. 
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through order a^. As the first step in this derivation, we rewrite this action using the covariant hypercubic fermion 
fields defined in Eq. ( p4| ) in which UA{y), the average of the link products along the shortest paths from site y to site 
y + A, is defined as: 



^^(y) = IT II U^,{y)^^U,{y+A^fL)^'^Up{y + A^fi+A,i))^^U^{y + A^^+A,i)+App)^% 



(A2) 



where the summation is on all permutations of {^vpa). We define the hypercubic first and second order covariant 
derivatives as 



1 



4a 



Df^^Aiy) - T- My)Ut.iy + P)VA{y + 2fi) - UUy - fi)Ul{y ~ 2fi)ipA{y ~ 2fi) 



(A3a) 



^^.VAiy) = -^ [Uf,{y)U^{y + fi)ipAiy + 2fi) + C/^(y - /i)f/^(y - 2fi)ipAiy - 2p) ~ 2^A{y) 



Then, with no approximation, we can rewrite Eq. (Al) as 



(A3b) 



SF = {2afY.Y.Y. VA{y)Ul,®I) 

y AB fj. ^ 



1 
AB 2 



wlP{y;AB)+w(^Hy;AB) 



-(75(8^5^)^5 



+(7m ® I)ab 

-(75«)C5m)aB 



ir 

2 

a 

2 
a 



w^^\y;AB)-w^^Hy;AB) 



i''< 



w]r'iy;AB)-w\fHy;AB) 






D,. 



D,. 



A„ 



wl'\y;AB)+w^^Hy;AB) 



+ h^^I)AB^ 
1 
4a 



u;W {y; AB) + ^f^) (y. AB) - w^^^ (y; AB) - u;(f ) ^y. ^b) 



- (75 ® ^5,)ab J: ^i'^ (y; ^B) + w^^^ {y, AB) - w^'^^ {y- AB) - u;(f ^ (y; AB)\ | tpe {y) 

+(2a)'^5m ^^A{y)vA{y)- 

V A 



(A4) 



The two hypercubic vectors A and B satisfy the delta function 5{A+B+ji)^ and the "bar" means modulo 2. The four 
closed-loop link products are defined as: 



w';l\y-AB) - 5A,iUA{y)Up{y+A)ul{y+2fL)Ul{y+p)Ul{y), 



(A5a) 



w^;^\y-AB) ^ 5A,oUA{y)Ul{y-2fL+B)ul{y-2p)U^{y-2fL)U^{y-^x), 



(A5b) 



wf\y-AB)^5A,iUA{y)Ul{y + B)Ul{y), 



(A5c) 



w'^^\y-AB) = 5A,oUA{y)Up{y+A)ul{y). 



(A5d) 



If we represent UA{y) by a doubled line which starts at site y and ends at site y + A, and U\{y) by a doubled line 
which starts at site y + A and ends at site y, then some typical loops can be shown by four figures (Figure |l| to ^. 

To expand the links in powers of "a" , we take advantage of the parallel transporter from x to x + jl to define the 
gauge field Af^ix) by the path ordered exponential 



Ufi{x) — Pexp< igofl / dTAi_i{x + rfi) 



(A6) 



When a ^ 0, we have 
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D^, = D,, + -a^Dl + Oia''), (A7a) 



A,=Dl + Oia'). (A7b) 

where D^ = 9^ + igo^f^ is the continuum covariant derivative. Expanding the w's in powers of "a", we get 

i[^W+,«f]=l + 0(a3), (A8a) 



2' 



I H^ - <^] = ^^50^2 E ^''^''M + 0(a3), (A8b) 






^hW+.«(^)-«;f-<^)] 2^ 



)] = ii5oa' E A, [D^, F,^] + 0{a^). (ASd) 

Finally, we get the classical continuum limit of the staggered fermion action as Eq. 



APPENDIX B: CALCULATION OF THE TREE-LEVEL COEFFICIENTS OF FOUR-FERMION 

OPERATORS 

We evaluate the amplitude represented by the graph shown in Figure |5| for the case of vanishing external momenta, 
pf — 0. The amplitude is given by: 

C08[^{7:c+7:D)MC+D + M+f}^). (Bl) 

where a, 6, c, and d are color indices, and fy^ is a hypercubic vector whose i^'s component is 1 only ii i> < ii. ttm is the 
momentum propagated by the gluon and the quantity tt^.^ is defined as 

7^M = ^y2sm^(M,7T/2). (B2) 

The primed summation on M is for all hypercubic vectors M with M^ = 0. 
After some algebra and using the notation of ref ||lj] , we get 



2 






^T.^^^)abWI7}cd}: (B3) 



with 



(7s ® Cf)^s = E l(-l)^''' i^s ® ^f)cd W^)""-^- (B4) 



CD 
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Then the action differs from the continuum by a four fermion term: 

^^^l H ll^AB'hDr{^A)x\^B)r{^c)X%^D). (B5) 



2 

ABCD abed 



Written in terms of the hypercubic fields, we get 



AS = -|(^^i2 + y^i4 + ^^13). (B6) 



The counterterm is the opposite of the above term. Hence we get the coefficients hsted in Eq. (pG 
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y 



y+2/i 



FIG. 1. A typical graph of w]^ {y;AB). It starts from y, through y+A, y + 2fi + B, j/+2/i, y+fJ., and then comes back to y. 
The two hypercubic vectors A and B satisfy the condition A^ = 1, Bf^ — 0. 



y-2/i+B 




y-2/i 



y 



FIG. 2. A typical graph of wj^ {y; AB). It starts from y, through y + A, y — 2fi + B, y — 2/i, and then comes back to y. The 
two hypercubic vectors A and B satisfy the condition A^ — 0, B^ = 1. 




y+A 



FIG. 3. A typical graph of w}, {y; AB). It starts from y, through y+A, y+B, and then comes back to y. The two hypercubic 
vectors A and B satisfy the condition A^ — 1, B^ — 0. 
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FIG. 4. A typical graph of wj, (j/; AB). It starts from y, through y+A, y+B, and then comes back to y. The two hypercubic 
vectors A and B satisfy the condition Afj, — 0, B^ — 1. 




d 



FIG. 5. The Feynman graph which generates four-fermion operators at the tree-level. 
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